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ABSTRACT 

Light curve observations of the double-mode Cepheid TD Cas 
obtained by 10 different sets of observers on several photometric 
systems over a tlnffi span of 67 years have been carefully studied to 
determine the fundamental and first overtone periods and their ampli- 
tudes on the V magnitude scale. The presence of a second overtone 
radial pulsation is discussed, and it is concluded that a previous 
detection of this mode was spurious due to the lack of a proper zero 
point correction for two groups of observations. The amplitudes of 
the two modes are shown to possibly vary during the entire observing 
period with the fundamental mode amplitude of 0.69 ^ 0.03 and the 
overtone amplitude decreasing about 0.2 or 0.3 magnitude. If this 
Cepheid displays the two pulsation inodes because it is mode switching, 
this switching time scale might be less than a hundred years. 
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Our motivation for this investigation is the reported 
existence by Faulkner (1977) of a third radial pulsation mode 
in the double-mode Cepheid TU Cas. At least one other Cepheid- 
type variable, AC And,* and many of the 5 Scuti variables have three 
radial pulsation modes, so it is perhaps not surprising to find an- 
other. There are two problems, however, TU Cas appears to be like 
the dozen or so double— mode Cephelds with a first overtone to fxinda— 
mental mode period ratio £= 0.71, unlike AC And with a ratio of almost 
0.74. In addition the reported third period is abnormally long with 
respect to the primary and secondary periods, making it very difficult 
to explain in terms of reasonable stellar models. The AC And model 
is within the normal range of composition and stiructure according to 
Cox, King, and Hodson (1978). 

Another motivation for studying TU Cas is to see if the ampli- 
tudes of the two modes are changing with respect to each other over 
the 67 years spanned by the available observations. Double-mode 
Cephelds, which might comprise one- third of all short period Cephelds 
(< 5 days) according to Stobie (1977) , have yet to be adequately 
explained. Stellingwerf (1975) suggested that double-mode behavior 
resulted at temperatures near the red edge of the Instability strip 
where the two modes switched toward each other. Unfortunately, Hodson 
and Cox (1976) found no double-mode behavior at these cooler tempera- 
tures. Currently, the only other cause of double-mode Cephelds 
predicted by theory is mode switching at the transition line between 
fundamental mode pulsation to the red and first overtone pulsation 
to the blue. 


26 



Observational evidence tends to support double-mode Cepheid pvilsa- 
tion near the transition line. Results from Rogers and Gingold (1973) 
on U Tr A, Schmidt (1972) on TU Gas, Pel and Lub (1978), and Cogan 
(1978) all arrived at temperatures placing most of these variables 
near the transition line. However, if double-mode behavior is the 
result of mode switching, we would expect to find a complete range 
of amplitude ratios between the two modes. As the evidence currently 
stands, only AX Vel has a first overtone amplitude larger than the 
fundamental. 

Amplitude changes over a span of 67 years would support the 
rapid mode switching rates predicted by Stelllngwerf (1975) and by 
unpublished results at Los Alamos, whereas no change would be more 
consistent with a possible slow change or just a stable permanent 
mixed mode. 

The observations used in our TU Gas analysis are given in 

Table 1. The set used by Faulkner contains 290 points and consist 

* 

of the four groups ranging from 1946 to 1959 (denoted by ). We 
attempted a period search on 302 white light observations by Osterhoff 
(1957), but were imsuccessful due to poor phase distribution of this 
data. For the amplitude analysis, we added 111 observations starting 
from 1962 to the present, which Include 60 new visual magnitudes by 
Schmidt observed from 1976-78. For a more sensitive mode switching 
search, we converted observations by van Biesbroeck and Gasteels 
(1914) from visual estimates to the modem V magnitude system with 
reasonably good results. 
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TABLE 1 


0BSERVATIC3NS OF TO CASSIOPEIAE 


REFERENCE 

EPOCHS 

(JD 2400 0 004-) 

YEAR 

NUMBER OF 
MEASURES 

OBSERVING 

SYSTEM 

Van Blesbroeck & 

19198-19481 

1911-1912 

192 

Visual 

Casceela (191A) 





it 

Gordon & Kron 

32039-32163 

1946 

28 


(1947) 




5000 

Oosterhoff 

32424-32471 

1947 

265 

Am , . - . , 

(1957) 

32606-32615 

1948 

37 

white light 


32272-33608 

1950 

350 


^ * 
Worley & Eggen 

33146-33206 

1949 

11 

V , (P-V) 

(1957) 

34334-34369 

1952 

4 

e* E 


35063-35113 

1954 

102 means 



35388-35404 

1955 

74 means 


* 

Oosterhoff 

36751-36849 

1959 

45 

UBV ' 

(1960) 





Weaver, Stelmetz & 

36742-36858 

1959 

26 

UBV 

Mitchell (I960)* 





Wllllama (1966) 

37936-37940 

1962 

6 

UBV 

Kvee & Braun 

38219-38235 

1963 

14 

UBV 

(1967) 





Takase (1969) 

38776-39079 

1965 

20 

UBV 


39358-39433 

1966 

11 

UBV 

Schmidt 

43015-43126 

1976 

26 

UBVRI 

(private communication) 

43403-43500 

1977-1978 

34 


* 

Observations used by 

Faulkner (1977) 

• 



For our periodicity search we 

fitted the 

Faulkner 

data group 


by least squares to a sinusoidal series 

m 

s 1 + cos2‘ir(if t-<t>^) (1) 


made up of a single frequency plus m-1 harmonics taken over a speci- 
fied frequency range. The quality of the fit and, consequently, the 
significance of any frequency is measured by a, the standard deviation 

of the observations from the fit. When m » 1, the amplitude is related 
2 2 

to a by A =2 (Aa ) . Equation 1 can be expanded to Include two or more 
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freqviencles , not harmonically related, and their cross coupling terms. 
The luminosity variations for the triple mode case are fit by 


im 

<L> 


1 + 


m 



|l| + U|+|k|-l 


^ijk cos2:r[(ifQ + + kf2>t - 


( 2 ) 


The appropriate order of a fit is determined when the variance of 
fit, order m + 1, is not appreciably smaller than the variance of m*^^ 
order fit. In addition, the fit will yield small changes in amplitude 
and phase angle with increasing order as suggested by 

Fitch. From Equation 2 we can directly obtain the amplitude and shape 
of each natural frequency. 

Third order Fouriergrams on the Faulkner data group yield the 
same primary and secondary period values as Faulkner obtained. We 
used these values to construct the double-mode fit, and then searched 
the residual for any remaining periodicity. In Figure 1 we see, in 
the upper curve, the first order Fourlergram of the fourth order resi- 
diial which is identical to Faulkner’s results. The change in o of the 
central minimum corresponds to an amplitude of ~ 0.014 for II 2 ~ 1.25247 
days. The sldelobes are due to the gap of ~ 1450 days between the 1959 
and 1954-55 epochs. The lower curve is the Fourlergram of the fifth 
order residual, and, although the central minimum still persist, its 
depth gives an amplitude of ~ 0.009, down by 35% from the fourth order 
amplitude. • 

In Figure 2, we see the same Fourlergram over a larger frequency 
range, expanded in the direction of slightly higher, but more theoreti- 
cally plausible frequencies. In the fifth order residual, we find at 
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TU CAS POURIERCRAM OF 4TH.5TH ORDER RESIDUAL <2S0) 



FREQUENCY (C/OAY) 

Fig* 1* First order (single frequency) standard 
deviation of the fourth and fifth order 
residuals versus frequency using the basic 
290 data points. 

TU CAS FOURIERGRAM OF 4TH.5TH ORDER RESIDUAL <290) 



Fig. 2. First order (single frequency) standard deviation of 
the fourth and fifth order residuals versus frequency. 
This figure shows the same data as in Fig. 1 but over 
a larger frequency range. 
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least one frequency at ~ 0.816 ^ with a deeper minimum than Faulkner’s 
in the far left portion of the curve, suggesting that perhaps his longer 
period may be an artifact of data noise. 

To test the quality of the data, we calculated the standard devi- 
ation of each epoch, grouped by year, about a fourth-order double-mode 
fit through all 290 points. The results are given in Table 2. From 
the a’s in the first column, we see that the 1949 data have rather large 
scatter. The 1952 data may also be noisy, but with only four observa- 
tions, its sigma value may be spurious. But the 1949 data remain suspect, 
especially considering the small scatter of the 1954-55 Worley and Eggen 
data. This epoch cotxld be intrinsically noisy, or it may not share 
the same zero point with the other epochs. 

We were able to confirm Faulkner’s magnitude zero point adjustments 
between observers, and ran a subsequent check on the zero point between 
different observing epochs. Our check centered around Worley and Eggen 


TABLE 2 

EPOCH STANDARD DEVIATIONS 
FOURTH ORDER DOUBLE-MODE FIT 


EPOCH 

Ho Shift 
(290) 

1955 Shift 
(290) 

1949 Shift’ 
(290) 

1955 + 1949 Shift 
(290) 

No Shift - 1949 
(279) 

1955 Shift - 1949 
(279) 

ALL (290) . 

0.0283 

0.0271 

0 .0267 

0.0254 

0.0231 

0.0214 

GK 1946 (28) 

.0260 

.0246 

.0257 

.0234 

.0252 

.0234 

WE 1949 (11) 

• 0835 

.0841 

.0713 

.0724 

.0000 

.0000 

WE 1952 (4) 

.0567 

.0553 

.0544 

.0552 

.0552 

.0555 

WE 1954 (102) 

.0205 

.0180 

.0212 

.0187 

.0208 

.0180 

WE 1955 (74) 

.0156 

.0095 

.0151 

.0089 

.0145 

.0077 

0+WSM 1959 (71) 

0.0323 

0.0332 

0.0308 

0.0314 

0.0299 

0.0307 


1955 Shift AV - -0.044 


1949 Shift AV - 0.066 
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data, as It is the only group with more than one epoch. We performed 
various trial shifts on the three epochs In this group In an attempt 
to reduce the relative and total standard deviation. 

From the third column, we see that by shifting the 1955 epoch by 
AV " - 0.044, the total o is reduced, and a for this epoch is down 
almost by a factor of two. A AV change of + 0.066 in the 1949 data 
also gave some reduction in the total sigma, and in the fifth column, 
we see that the combined effects of shifting both 1949 and 1955 epochs 
results in an ever smaller total deviation. 

Removing the 1949 epoch from the unshifted 290 points produces , 
by itself, a larger reduction in total sigma than does any of pre- 
vious shifting. 18% reduction in O for a 4% point reduction) . 
Removing this epoch and shifting the 1955 epoch gives the largest 
reduction in sigma, almost 24% lower than for the original Faulkner 
data. 

In Figure 3 we have Fouriergrams of fourth order fit residuals 
in the region of Faulkner's third period for each trial listed on the 
table. The uppermost curve is the same Fouriergram as you saw in 
Figure 1, showing the central minimum surrrounded by aliases. By 
shifting 1955 epoch this pattern almost disappears (second curve from 
top). However, in the next curve, we see that 1949 epoch shift alone, 
produces only small changes from the original pattern of the top curve. 
The fourth curve shows the combined effect of the two shifts, and is 
similar to the second curve because of the large effect of the 1955 
epoch shift. The fifth curve results from removing 1949 epoch from 
the data, and in spite of a substantial reduction in O’, the original 
pattern still persists. This pattern again disappears when we 
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TU CAS FOURIERCRAM OF 4TH ORDER RESIDUALS 



the several fourth order fit residual cases of 
Table 2 versus frequency. 


introduce the 1955 data shift as seen in the bottom curve. 

Expanding the frequency range in Figure 4, as we did before, 
shows that shifting 1955 data removes the cyclic pattern of curves 
one, three and five, which corresponds to the 291 day interval 
between the 1954 and 1955 epochs. Faulkner's proposed third period, 
at the far left, appears to be only a deeper meciber of this pattern 
that disappears with the appropriate zero point shift. 

Fitch noted the amplitude given by Faulkner's limited triple 
mode fit was a factor of two larger than the amplitude predicted by 
the depth of the minimum in the Fouriergram (about 0.025 to 0.014). 
We can explain this discrepancy in Figure 5 with two plots of funda- 
mental and first overtone prewhitened data and mean curves as a 
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TU CAS FOURIERCRAM OF 4TH ORDER RESIDUALS 



FREQUENCY (C/OAY) 

Fig. 4. First order (single frequency) standard deviation of the 
several fourth order fit residual cases of Table 2 versus 
frequency. This figure shows the same data as in Fig. 3 
but over a larger frequency range. 

function of phase. In both cases we have tised the original, uncorrected 
Faulkner data set. We produced the upper plot by removing all the 
triple-mode fit terms except the II 2 sine term, and phasing the resi- 
dual to il 2 . This plot is identical to Faulkner's (Figure 4, 1977) and 
confirms his amplitude for II 2 given by the triple— mode fit. The lower 
plot is the phased residual of the fourth order-double mode fit, with 
the amplitude estimated by a spline fit. It gives an amplitude value 
of ~ 0.017 which is much closer to the predicted value of 0.014, with 
the difference of 20% due to uncertianties in the spline fit such as 
the number of knots and order. In general, we have noticed when any 
new natural mode (real or spurious) is introduced into the multimode 
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TU CAS UPP£R-LTMF-AMPL=0.02S/L0W£R-0MP-AMPL=0.017 



the points obtained by prewhitening by all the lU 
and ITi terms of his fit- The lower panel gives the 
variation of the residuals of a double-mode fit with 
the proposed phase together with a spline fit 
through the points. 

fit, it interacts in varying degrees with the other natural modes, 
resulting in a substantial amplitude for this mode. Hence, one 
cannot reliably determine the significance of a variation by its 
amplitude size given by multimode fitting. 

In searching for amplitude changes in TU Gas, we selected the 
longest time base possible with the hope of revealing the more 
conspicuous variations in modal amplitude. A summary of our results 
is given in Table 3. 

The 1911-12 epoch by van Biesbroeck and Casteels was converted 
from a visual estimates system to V magnitudes. Initially, our 


35 




TABLE 3 


MODE AMPLITUDES OF TU CASSIOPIAE* 


Epoch 

No, of 

Observations 

^0 


0 

Aj^/A| 

1911-12^ 

192 

.70 

.45 

.13 

.64 

1911-12^''’ 

192 

.72 

.54 

,19 

.74 

fjf 

1911-12 

192 

.86 

.53 

.14 

.62 

1946-1977 

390 

.71 

.27 

.04 

CO 

19A6-1959 

279 

.71 

CM 

.03 

.41 

1962-1977 

111 

.67 

.24 

.03 

.35 


Third order double-mode Fourier Series fit 

-U 

Adopted magnitudes by Van Biesbroeck and Casteels 

Revised magnitudes from HD catalog 
tft UBV magnitudes by Henden (private communication, 1978) 


primary uncertainty in converting these observations resulted from 
the absence of any recorded UBV measurements of the six comparison 
stars used for this study. Van Biesbroeck and Casteels estimated 
these magnitudes from a relative scale used to determine the variable's 
changing magnitude. We attempted to iiqjrove on these magnitudes using 
the best revised estimates from the HD catalog (as suggested by 
Bidelman, private communication). Eventually we were able to obtain 
UBV magnitudes for the six comparison stars from Ame Henden at 
Indiana University (private communication, 1978) 

The ratio of first overtone to fundamental mode amplitude 
(A^/Aq) for this earliest epoch averages ~ 0.67 and is appreciably 
higher than for the 1946-1977 era (~ 0.38). To check this last ratio 
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we split the modern era into two ~ 15 year intervals to determine 
if this trend still persisted. The decrease of the ratio from 0.41 
to 0.35 is suggestive, but probably within the range of uncertianty. 

Due to the relatively larger uncertainty of the 1911-12 data, 
we are unable to conclude that the first overtone mode really decayed 
substantially since 1911. However, observations after 1945 span 
nearly 32 years and comprise a long time base of quality UBV measures. 
With the inclusion of new observations the decrease of within 

the modem epoch may possibly be confirmed. 

Our conclusions are that TU Gas is not a triple-mode pulsator, 
and that maybe the first overtone pulsations have decayed since 1911. 

We wish to thank R. F. Stellingwerf for locating the van Biesbroeck 
and Casteels data and providing it to us. 
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Discussion 


Wesselink ; How do you identify a mode? How do you know whether it is the 
fundamental or an overtone? 

Hodson ; Hopefully, the one with the longest period with largest amplitude 
is the fundamental. Sometimes that turns out to be the 1st or 2nd overtone. 

But that's the nomenclature I use. 

J. Cox ; In 1977, I think, Stobie pointed out that there was no observational 
evidence for amplitude changes over at least a 50-year time span. Do you 
think this is because no one has looked, or is there more substantial evidence? 

Hodson ; The evidence we have suggests that this might be happening. We have 
32 years of observations in the "modern" era. We need quite a few more 
"modern" observations to resolve the problem. 
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